Giant Pulses from the Millisecond Pulsar B1821— 24 
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ABSTRACT 

We have carried out a survey for 'giant pulses' in 5 millisecond pulsars. We detect individual 
pulses from the high E pulsar PSR B1821— 24 with energies exceeding 50x the mean pulse energy. 
These giant pulses are concentrated in a narrow phase window coincident with the power-law 
non-thermal pulse seen in hard X-rays. This is the third example of the giant pulse phenomenon. 
It supports the idea that large B fields in the outer magnetosphere are critical to the formation 
of such pulses and further suggests a direct connection between giant pulses and high energy 
emission. 



Subject headings: pulsars: individual (PSR B1821— 24) 
1. Introduction 

Despite three decades of intensive study, the 
mechanism producing pulsar radio emission is 
poorly understood. Fluctuations in the inten- 
sity of the radio radiation provide important con- 
straints on plausible mechanisms. Single-pulse 
studies of bright pulsars detect a variety of pat- 
terns in the intrinsic intensity fluctuations, in- 
cluding nulling and drifting phenomena. The dis- 
tribution of integrated pulse energies, however, 
has only a modest dispersion, typically forming a 
Gaussian of width < (E) , the mean pulse energy 
(eg. Manchester & Taylor 1974). In contrast, the 
Crab pulsar PSR B0531+24 emits pulses with en- 
ergies > 20 x (E), extending up to > 2 x 10 3 (E) 
(Lundgren et al. 1995) which were instrumental 
in the original detection of the Crab (Staclin 
& Reifenstein 1968). These Giant Pulses (GP) 
are typically broadband (Moffett 1997, Sallmen 
et al. 1999) and of short duration, with widths of 
order a few /is and structure down to 10ns (Han- 
kins 2000). They are localized to the main and 
interpulse phase windows and follow an intensity 
distribution best characterized as a power law with 
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index of ~ 3 — 3.5. 

The discovery of similar pulses from the mil- 
lisecond pulsar PSR B1937+214 (Sallmen & 
Backer 1995, Cognard et al. 1996) was surpris- 
ing. The pulses are extremely short (r < 0.3/is) 
events confined to small phase windows trailing 
the main and inter- pulses, again with an ap- 
proximately power-law distribution of pulse ener- 
gies (Kinkhabwala & Thorsett 2000). Since PSR 
B1937+214 is the only known radio pulsar having 
an estimated magnetic field at the light cylinder 
larger than that of the Crab, it has been suggested 
that this is a key parameter controlling the giant 
pulse phenomenon (Cognard et al. 1996). 

Johnston et al. (2001) have recently found that 
the Vela pulsar produces pulses with a very wide 
distribution of peak fluxes. It is not clear if these 
events are related to true giant pulses; the largest 
such pulses observed to date have E < 10(E), 
but these narrow pulses have peak fluxes exceed- 
ing 40 x the integrated peak intensity. Relatively 
few individual pulses have been examined, and the 
extended tail in the E distribution of these 'giant 
/x-pulses' may continue into the true GP regime. 
The pulses are localized to a phase window off of 
the integrated pulse emission, are of short dura- 
tion and are highly polarized. 

To explore the connection between the Crab 
and PSR B1937+214 GPs and the large individ- 
ual pulses of 1 — 3 x 10 4 yr Vela-like pulsars, we 



obtained fast time sample data for several young 
and millisecond pulsars. We report here evidence 
for giant pulses in the MSP with the next high- 
est known light cylinder field, PSR B 182 1-24, 
along with limits on GP emission in other high- 
field MSPs. In a companion paper (S. Johnston 
& R.W. Romani 2001, in preparation) we discuss 
large individual pulse energies in young high field 
pulsars. 

2. Observations 

Observations were made on 20-22 May 2001, 
with the Parkes 64-m radio telescope. We used the 
center beam of the 21-cm multi-beam system at an 
observing frequency of 1517.75 MHz. The receiver 
has a system equivalent flux density of 30 Jy on 
cold sky. The back-end consisted of a filtcrbank 
containing 512 channels per polarization each of 
width 0.5 MHz for a total bandwidth of 256 MHz. 
The polarization pairs are summed, each output 
is then sampled at 80 /is, one-bit digitized, and 
written to DLT for off-line analysis. 

The data were then de-dispersed at the pulsar's 
nominal dispersion measure. The mean and rms 
of groups of 8192 samples (0.65 s) were examined 
and those which showed obvious signs of inter- 
ference were discarded. The data could then be 
folded synchronously with the pulsar's topocentric 
period to produce a pulse profile. Our nominal 5- 
cr sensitivity in 80 fis is 1.3 Jy. In practice, even 
after clipping, we experienced substantially larger 
background fluctuations and the rms exceeded our 
expected rms by a factor of 1.5-3. 

To confirm that we could detect conventional 
GPs we made short observations of the Crab pul- 
sar, during which GPs were detected with high 
significance. Short integrations with the Vela pul- 
sar off-axis, producing an effective continuum flux 
of only 0.6 mJy, confirmed that we could detect 
the largest individual Vela-type giant ^-pulses at 
these faint continuum flux levels. Finally, and 
most relevant to the present study we observed 
PSR B1937+214, obtaining 1784 s of integration 
(1.15 x 10 6 pulses). Our sampling provides only 
19 bins across the pulse profile, but we clearly de- 
tect the GP distributions in both the main and 
intcrpulsc components. In both components the 
giant pulse peaks occur approximately one bin af- 
ter the corresponding peak of the integrated pulse 



profile. The largest main component GP obtained 
had an energy 445 Jy/its, the third largest had 230 
Jy/is. The phasing and intensity distribution are 
well matched to the 1.4 GHz results reported by 
Kinkhabwala & Thorsctt (2000). 
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Fig. 1.— Giant pulses from PSR B1821-24. Up- 
per panel: the peak flux of the brightest pulses 
with adjacent bins > 4cr, points at the peak phase. 
Our integrated profile (arbitrary scale) is overplot- 
ted for reference; 1.2 periods are shown. Lower 
panel: Integrated energy of the brightest GP. We 
label the components on the reference pulse pro- 
file and mark (at top) the 4-bin windows used to 
measure component energies. 



3. Results 

PSR B1821-24 has a period of 3.05 ms and a 
dispersion measure (DM) of 119.8 cm~ 3 pc. This 
DM leads to a smearing of 142 /is across our 0.5 
MHz filter bands, limiting our resolution and en- 
suring that any giant pulse from the source will 
be smeared across 2-3 time samples. We obtained 
3 x 3600 s of data on this pulsar starting at MJD 
52049.7; after excision of 0.65 s averages showing 
strong interference, we retained 9943 s of data or 
«3.3xl0 6 pulses. Because of the expected DM 
smearing, we first searched the de-dispersed data 
stream by tagging pulses with 2 or more consec- 
utive samples exceeding 4o\ Our data contained 



16 such pulses. The brightest of these, both in 
terms of peak sample intensity and integrated en- 
ergy were confined to a narrow window in the first 
component of the integrated pulse profile (Figure 
1 , which shows the brightest 8 pulses by each mea- 
sure). During these observations PSR B1821— 24 
had a continuum flux of «3 mJy and an average 
pulse energy of 9.3Jy/zs. The largest individual 
pulse observed had E « 755Jy^s which is 81(E). 
This pulsar has a complex pulse morphology. 
We adopt the naming convention of Backer & Sail- 
men (1997), in which PI has a steep spectrum, 
dominating the pulse at low radio frequencies and 
the second narrow pulse P2 is stronger at short 
wavelength. The broad trailing component P3 
seems to have a similar spectrum to PI. Backer 
& Sallmcn report that at 1.4 GHz the pulsar oc- 
casionally spends many hours in a mode with P2 
decreased relative to the other components by a 
factor of 4. We do not see such dramatic varia- 
tions, but confirm that P2 is variable, decreasing 
in flux by 25% relative to PI and P3 over our first 
hour of integration and stabilizing afterwards. 
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Fig. 2. — Above: the average profile of the 16 
GP candidates (full line histogram) compared with 
the average profile of all pulses, magnified by 40 x 
(dashed line curve). Below: the RXTE pulse pro- 
file with the relative phase determined by Rots 
ct al. (1998). 

The pulsar is also detected in the X-ray, again 



with a complex pulse profile. Careful timing (Rots 
et al. 1998) using the Rossi X-ray Tming Explorer 
(RXTE) establishes that the narrow, hard spec- 
trum X-ray component is associated with PI, but 
lags significantly by 60 ± 20/is. A second broad, 
and apparently softer, X-ray peak lags P3 by 
250±60^s. Figure 2 compares our integrated pulse 
profile and the profile obtained by averaging all 16 
pulses with adjacent samples > 4c. These pulses 
have been re-sampled to common phase bins, with 
38 points across the light curve; each bin corre- 
sponds to about one original sample. The 'GP' 
lightcurvc PI emission lags that of the integrated 
pulse profile by approximately 80/xs (one bin). No 
other significant GP feature is seen. For compari- 
son we plot the RXTE light curve of Rots et al. , 
with the relative phasing determined in that study. 
DM smearing and re-sampling to fixed phase win- 
dows can fully account for our GP pulse width. 
Thus our GP emission appears to be unresolved 
in phase and coincident with the hard power law 
high energy pulse. 

We can obtain some further information on the 
GP distribution by measuring pulse intensity dis- 
tributions in windows of pulse phase. For each 
pulse we defined a 10 sample off pulse window (<p = 
0.63-0.85 in Figure 1). The individual components 
were estimated by the three 4-sample windows in- 
dicated in Figure 1 (Pl=-0.03 - 0.05; P2=0.24 
- 0.32; P3= 0.45 - 0.53). For each pulse we de- 
termined the background level from the median of 
the samples in off pulse region of the current and 
preceding pulses. The individual component en- 
ergy was measured by an integration of the energy 
in the appropriate four sample window. The cu- 
mulative histograms of the three components arc 
shown in Figure 3. Component PI shows an ap- 
preciable excess of high energy pulses, extending 
somewhat below our selected giant pulse thresh- 
old. Note that due to DM smearing and our lim- 
ited resolution, some pulse energy will be lost from 
these 4-sample windows, so these E values under- 
estimate the full pulse energy 

Since only our brightest GPs are well separated 
from the underlying pulse and noise fluctuations, 
we can put little constraint on their intensity dis- 
tribution. Power law indices of w 3 — 5 ade- 
quately describe the GP excess. Indeed an addi- 
tion of an E~ 5 power law to the P3 flux distribu- 
tion nominally matches the PI distribution quite 
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Fig. 3. — Cumulative probability distributions of 
pulse energy in our three selected pulse component 
windows. The probability histograms are contin- 
ued as smooth lines at low pulse energies. There 
P2 is the highest, as expected from the average 
pulse profile. At large energies there is a clear 
excess in component PI (full line histogram). A 
simple power law for the excess GP distribution is 
plotted with the dashed line. 



well (dashed smooth line). As for the 1.4 GHz 
PSR B1937+214 data of Kinkhabwala & Thorsett 
(2000), it is likely that noise fluctuations and nor- 
mal pulses steepen the estimated power law index 
this close to the noise floor. Longer observations 
with greater sensitivity will be needed to probe the 
true GP intensity distribution. We expect about 
1/h (P « 8.5 x 10~ 7 ) with E > 28(E). This 1.4 
GHz GP rate is, in fact, roughly twice that of 
PSR B1937+214, which has P (28(E)) « 4 x 10~ 7 . 



If our GP are as narrow as those of PSR 1937+214, 
we expect the brightest pulse in an hour to reach 
a peak flux of > 1300 Jy. 

4. Limits for other MSPs 

We also observed four other millisecond pul- 
sars with high inferred light cylinder magnetic 
fields, -Blc- In Table 1, in i?LC order, the pulsars 
that we observed are listed along with the mean 
pulse energy, number of pulses observed and the 
limit to which giant pulses were not seen. These 
four pulsars have -Blcs substantially lower than 
that of PSR B1821-24. For PSRs J0034-0534, 
J2129-5721 and J1911-1114, the low flux at our 
observation epoch precluded significant bounds on 
GP emission. The last entry of the table lists the 
factor by which a pulse must exceed (E) to be 
detected in our study. Indeed, if the cumulative 
intensity distribution followed Prob(> E) oc E~ 3 
and the GP rate were the same as for B1821— 24, 
we would have required 10 5 , 58 and 18 million 
pulses, respectively, to see a giant pulse in these 
three pulsars at our achieved threshold. Consid- 
ering further the likely decreased GP frequency at 
lower Blc i ° ur limits do not significantly constrain 
GP emission in these sources. On the other hand, 
for J0613— 0200 (with a similar period and flux 
density to PSR B1821-24 but B LC lower by 5), 
our limits do begin to probe B 182 1—24 -like GPs. 
For the GP distribution assumed above with the 
same probability as in B1821— 24, we would have 
expected about 15 GP events. In fact, the largest 
few pulses do occur in 0.1 of phase ahead of the 
integrated pulse profile, suggesting that more sen- 
sitive observations might uncover a GP distribu- 
tion in this source at < 1/10 the frequency of the 
events seen in B1821— 24. 



Table 1. MSP Giant Pulse Candidates 
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5. The Giant Pulse Phenomenon and 
High-Energy Emission 

The addition of B1821-24 to the GP family 
supports the idea that large magnetic fields in the 
outer magnctosphere are required for this behav- 
ior. However this object points to an additional 
common feature: as for the Crab pulsar, the GPs 
appear isolated to narrow windows of pulse phase 
showing high energy pulses, with hard power law 
spectra. In the case of the Crab these X- and 7-ray 
pulses are believed to arise in an outer magncto- 
sphere acceleration gap along the boundary of the 
open zone (Cheng, Ho & Ruderman 1986; Romani 
& Yadigaroglu 1995). These Crab-type outer gaps 
are believed to have 7 — 7 pair production main- 
tained by a dense bath of target soft synchrotron 
photons from the gap itself. High B fields (espe- 
cially near the null charge surface) should enhance 
this synchrotron emissivity, and hence the rate of 
secondary e ± production. This dense pair plasma 
produces the narrow hard X-ray pulses, and we 
can speculate that the high densities promote the 
instabilities that create enhancements in the par- 
ticle coherence and hence the giant radio pulses. 

This picture provides a few predictions. PSR 
B1937+214 shows narrow X-ray pulses in ASCA 
data (Takahashi et al. 2001), with a phase close 
to the radio interpulse. We expect that a high 
resolution X-ray profile with good absolute radio 
phase comparison will show these X-ray pulses to 
align with the GP phase, lagging the integrated 
radio profile pulse peak. It is not clear why the 
interpulse should be stronger at high energies than 
the main pulse. An additional strong GP candi- 
date not accessible to our survey is the millisec- 
ond pulsar PSR J0218+4232 (5 LC = 3.1 x 10 5 
G). Comparison with the high energy pulse pro- 
file (Kuipcr et al. 2000) suggests that GP emission 
may be found in the second radio pulse, and pos- 
sibly in the > 100 McV pulse phase preceding the 
first radio peak. 

The detailed dependence of the GP frequency 
on Blc is not clear. If the critical quantity is the 
field near the outer gap base at the null charge 
surface, where pair production is expected to be 
most intense, then magnetic inclination a can have 
a strong effect with B w i?Lc[l-5tan(a)] 6 for a 
static dipole field and a > 30°. Consider the 
case of the Crab pulsar, which has a frequency 



of E > 20(E) pulses roughly 10 4 x larger than 
that of PSR B1937+214 or PSR B1821-24 at 1.4 
GHz. Polarization modeling (Romani and Yadi- 
garoglu 1995) and the recent CXO wind torus im- 
ages suggest that the Crab has a « 80°. Baker 
and Sallmen (1997) found a plausible polarization 
sweep fit of a w 50° for PSR B1821-24. For this 
angle the nominal field at the base of the gap is 
suggestively 10 4 x lower than that of the Crab. Of 
course, the magnetosphere of a MSP is likely far 
from vacuum dipole, but this does underline that 
a can have a strong effect on the ordering of GP 
candidates by B. 

A better criterion for the presence of strong 
GP activity may be the presence of a narrow 
hard spectrum X-ray pulse, indicating a dense 
pair plasma in an active outer gap. Most known, 
young 7-ray pulsars are Vela-like with pair pro- 
duction controlled by thermal surface fluxes (e.g. 
Romani 1996). These do have non-thermal X-ray 
pulse components, but these are much broader and 
the emission represents a much smaller fraction of 
the spindown luminosity. GP distributions might 
therefore be much weaker in these pulsars. The 
occurrence of GP in other young Crab-like pul- 
sars and the connection with the large pulses from 
millisecond pulsars described here can only be ad- 
dressed by further sensitive single pulse studies. 

The Parkes telescope is funded by the Com- 
monwealth of Australia for operation as a National 
Facility managed by CSIRO. The 512 channel fil- 
terbank system used in the observations was de- 
signed and built at the Jodrell Bank Observatory, 
University of Manchester. Observing software was 
provided by the Pulsar Multibeam Survey group. 
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